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In the framework of the instanton vacuum model we make expansion over the current mass m 
and number of colors Nc and evaluate 0(1/Nc, m, m/Nc, m In m/A^c)-corrections to the dynamical 
quark mass M, the quark condensate {qq), the pion mass and decay constant F^. There are 
several sources of these corrections: meson loops, finite size of the instanton distribution and the 
quark-quark "tensor" interaction terms. In contrast to the expectations, we found that numeri- 
cally the 1/A^c-corrections to dynamical mass are large and mostly come from meson loops. As a 
consequence, we have large 1/A^c-corrections to all the other quantities. To provide the values of 
Fn{m = 0), {qq{m = 0)) in agreement with xPT, we offer a new set of parameters p, R. Finally, we 
find the low-energy SU{2)f chiral lagrangian constants h,l4 in a rather good correspondence with 
the phenomenology. 
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I. INTRODUCTION 

The spontaneous breaking of chiral symmetry (SxSB) 
is one of the most important phenomena of hadron 
physics. It defines the properties of all the light mesons 
and baryons. Using the general idea of chiral symme- 
tryjthe Chiral Perturbation Theory (xPT) was proposed 
in [l| for evaluation of the QCD hadronic correlators at 
low-energy region, where the expansion parameters are 
light quark current masses m and pion momenta p. The 
basic tool is the phenomenological effective lagrangian, 
which has a form of the infinite series in these parame- 
ters p and m. Naturally, the low-energy constants (LEG) 
of the series expansion are not fixed. Up to now they 
were extracted only from the experimental data. Recent 
progress in lattice calculations provide us with rough es- 
timates of LEG. The main problem of lattice evaluations 
are the still-large pion masses available on the finite 
size lattices. 

QGD instanton vacuum model, often referred as the in- 
stanton liquid model, provides a very natural nonpertur- 
bative explanation of the SxSB [3, Q. It provides a con- 
sistent framework for description of the pions and thus 
may be used for evaluation of the LEG. 

Quasiclassical considerations show that it is energet- 
ically favourable to have lumps of strong gluon fields 
(instantons) spread over 4-dimensional Euclidian space. 
Such fields do strongly modify the quark propagation due 
to the t'Hooft type quark-quark interactions in the back- 
ground of the instanton vacuum field. This background 
is assumed as a superposition of N^^ instantons and N_ 
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antiinstantons 



N- 



(1) 



where ^ = (p, z, U) are the (anti)instanton collective co- 
ordinates ~ size, position and color orientation. The most 
essential for the low-energy processes are the would-be 
quark zero modes, which result in a very strong attrac- 
tion in the channels with quantum numbers of vacuum, 
appearance of the quark condensate and generation of 
the dynamical quark mass (see the reviews 

The main parameters of the model are the average 
inter-instanton distance R and the average instanton size 
P- 



The estimates of these quantities are 
p ~ 0.33 /m, i? ~ 1 fm, (phenomenological) [1], 



(2) 



p ~ 0.35 fm, R ~ 0.95 fm, (variational) 0, 

p ~ 0.36 fm, R ~ 0.89 fm, (lattice) @, 0, i, i, [I^ 

and have ~ 10 — 15% uncertainty. 

Recent computer investigations |ll| of a current mass 
dependence of QGD observables within instanton liquid 
model show that the best correspondence with lattice 
QGD data is obtained for 



p~0.32/m, i?~0.76/m. 



(3) 
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While in the real world the number of colors is iVc = 3, 
since the pioneering work of t'Hooft [T^ it is assumed 
that one can consider A^c-counting as a useful tool, i.e. 
take the limit Nc oo and neglect the 1/A^c-corrections. 
In the instanton vacuum model iVc-counting is naturally 
incorporated. The phenomenological set ^ is popular 
since in the leading order (LO) on Nc it reproduces rea- 
sonable values for most of the physical quantities. This 
leads to rather consistent description of pions and nucle- 
ons in the chiral limit. 
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The main purpose of this paper is evaluation of 
0{1/Nc, m, m/Nc, m/Nc hi TO)-corrections to different 
physical observables. There are several sources of such 
NLO corrections: 

1. At pure gluonic sector of the instanton vacuum 
model the width of the instanton size distribution 
is 0{1/Nc)- The account of the finite width leads 
to rather small corrections [l^. In the following 
we will check the accuracy of i5-function type of 
the instanton size distribution by direct evaluation 
of the finite width corrections. 

2. The back-reaction of the light quark determinants 
to the instanton vacuum properties is formally 
controlled by Nj /Nc-factov It does not 
sizably change the distribution over + but 
radically change the distribution over Nj^. — 
Any m/ = leads to (5-function type of the distri- 
bution . In the following we take = iV_ . 

3. There are the quark-quark tensor interaction terms 
which are l/TVc-suppressed and thus are absent in 
the LO effective action. These terms correspond to 
nonplanar diagrams in old-fashioned diagrammatic 
technique. 

4. The contribution of meson quantum fluctuations 
(meson loops) has to be taken into account. 

In the first few sections we study the role of the me- 
son loops which give the dominant contribution. At the 
end we also estimate the contributions of finite width of 
instanton size distribution, and above-mentioned tensor 
interaction term. 

We consider parameters p, R as free within their 15% 
uncertainty and fix them from the requirement -F',r,m=o = 
88MeV, {qq)m=o = {255MeVf with account of NLO cor- 
rections, as it is requested by xPT [ij. We found the 
values 

p = 0.350 fm, R = 0.856/m, (4) 

in agreement with the above-given estimates (j2l3|) . 

Note that though the evaluation of the meson loop 
corrections in the instanton vacuum model is similar to 
the earlier meson loop evaluations [l^, [l^, [3 
NJL model, there are a few differences which should be 
mentioned: 

1. As it has been already mentioned, the meson 
loop corrections are not the only sources of l/N^- 
corrections in the instanton model. 

2. Due to nonlocal formfactors there is no need to 
introduce independent fermion and boson cutoffs 
A/,Ab. The natural cutoff scale for all the loops 
(including meson loops) is the inverse instanton size 



3. The quark coupling constant is defined through 
the saddle-point equation in the instanton model 
whereas it is a fixed external parameter in NJL. 

In Section ini using the above mentioned assumptions 
and zero- mode approximation (see below), we calculate 
the quark determinant in the gluon background ([T]) and 
flavour vector w^, axial- vector a^, scalar s and pseu- 
doscalar p external fields beyond chiral limit. Since 
the zero-mode approximation leads to a violation of the 
flavour local gauge invariance, we formally restore it by 
introducing the path-ordered exponents. The price of 
such restoration is the path-dependence. We show, how- 
ever, that the quantities considered in this paper do not 
depend on the choice of the path (see Sec. IX)) . 

After that, we average the quark determinant over in- 
stanton collective coordinates as it was described in [l^, 
[20l . [2l| . using smallness of packing parameter (^)^ ~ 
0.01. This leads to the partition function Zn[v, a, s,p, to], 
which is a generating functional of various quark correla- 
tors. The partition function is represented as a path in- 
tegral over the fermions, which are in the present formal- 
ism constituent quarks with t'Hooft-like 2A^/-legs nonlo- 
cal interaction term with the coupling constant A, defined 
dynamically from the appropriate equation. In the fol- 
lowing we consider only Nf = 2 case. 

The simplest and the most important quark correla- 
tor is the quark propagator, since all the other observ- 
ables are sensitive to its properties. In section IIIII we 
calculate the dynamical quark mass M . We found in the 
instanton vacuum model that contrary to the large-iVc 
expectations, the meson-loop corrections to M[m) are 
large (~ 50%) even in the chiral limit. These corrections 
are intimately related to the equation for the coupling A. 
Such corrections are absent in NJL-type models. 

All meson-loop corrections to physical observables can 
be splitted into two parts: "indirect" (which are domi- 
nant and come from the dynamical mass shift) and "di- 
rect" (which are relatively small and thus obey 1/Nc- 
counting rules.) 

In sections ITVllV Al we calculate the meson-loop cor- 
rections to the quark condensate (qq), pion mass 
and decay constant Fy^. From 0{m, 1/Nc) corrections to 
the latter quantities in section IVIIII we extract the low- 
energy parameters ^3, 14 of the phenomenological Gasser- 
Leutwyler SU{2)f lagrangian [l|]. These two constants 
are especially important for chiral perturbation theory 
[2^ . [23I since they relate the physical observables F^^ , M^^ 
with their parameters. In section HXl we check that the 
chiral log theorems [Ull^, UHl are satisfied. In section 1x1 
we demonstrate that our results are path independent. 
Appendix |^ contains explicit expressions of the nonlocal 
vertices of vector and axial currents. In Appendix iBl we 
give details of evaluation of finite width corrections, and 
in Appendix [C] we study the structure of the correlator 
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II. LIGHT QUARKS PARTITION FUNCTION 
IN EXTERNAL FIELDS 

The hght quark partition function Z[m,v,a, s,p] with 
quark mass m and in external vector w^, axial- vector a^, 
scalar s and pseudoscalar p fields is defined as 

a, s,p, m] (5) 
= / DAfj_e~^ det{iD + im + V + a-f5 + s + pjr,). 



The basic assumption of the instanton model is that one 
can evaluate the integral in quasiclassical approximation, 
making expansion around the classical vacuum ([T]). The 
first evaluation of the partition function ^ was per- 
formed in Ref. [H, in the absence of the external 
fields and in the chiral limit. The main purpose of the 
present paper is the extension of the result of [l^, to 
the case of nonzero quark mass m and external v,a, s,p 
fields. While inclusion of the spin-0 fields s,p is triv- 
ial, spin-1 fields U/^,a^ deserve special attention due to 
the nonhomogeneous transformation under flavour rota- 
tions. Following [l^, [3l we split the quark determi- 
nant into the low- and high-frequency parts according to 
Det = Dethigh • Detiow and concentrate on the evaluation 
of Detiow, which is responsible for the low-energy domain. 
The high-energy part Dethigh is responsible mainly for 
the perturbative coupling rcnormalization. 

The starting point of our consideration is the zero- 
mode approximation formulated in p^ . [l^ . [26| 



p + A, 



1 ^ l^»,o)(^»,ol 
P 



im 



(6) 



Here the zero-modes $i^o are also functions of the instan- 
ton collective coordinates £^i. This approximation is good 
for small values of m (chiral limit) and indicates that the 
main contribution to the quark propagator is due to the 
zero-modes. The extension of ^ beyond the chiral limit 
was proposed in our previous works [U, [l^, [H, [2^ as 



where 

Q = -{<^0^\pSop\^Qi). 



(7) 



(8) 



The advantage of this approximation is that it gives cor- 
rect projection of Si to the zero-modes: 



S^\^0^) - — I'i'o,), ($0.1^. = ($0,1 — , 

im ira 



(9) 



while the similar projection of Si given by Eq. ([6]) is valid 
only at TO ^ limit. 

Recently the quark determinant Dct(p + Ai + im) ~ 
exp(r/) was calculated in the field of a single instanton 
and beyond the chiral limit [13, [3l| . The F / can be rep- 
resented asT f = — In TO/9 — 2Ts, where we take the scale 
equal to 1/p. Here the first term is due-to the zero- mode. 
At smah mp the F^ = 0.145873 + 0.5{mpy \n{mp) - 



0.0580(mp}f . First, independent on to term was calcu- 
lated at |32l |. It is especially interesting for us the to- 
dependent terms in F^. In the region mp ~ 0.2 (corre- 
sponding to the strange current quark mass) the third 
term is of order 0.2% of the second one and is negligible. 
So, the main problem is to reproduce the second term in 
Fg. Simple calculations with using the Eq. ([7]) lead to 
the Fg = (0.5 In m/9 + 0.0767) (mp)^. Again, at the region 
mp < 0.2 we may neglect by [mp)"^ term. So, our approx- 
imation ([7]) leads to the main (mp)^ In mp-term in Fg, 



which coincide with exact calculations [3C 
contrast to the zero-mode approximation i 
to r, = 0. 

We use the notations 

^ - S ^ 



S 



P + A- 
So^ 



V- 
1 



p + A, + V + ; 



l3ll | in a sharp 
2]) , which leads 



(10) 



P + V + ' 



So 



1 



p 



(11) 



for the quark propagators in the instanton background 
(fT|) or individual instantons Ai and in the external flavour 
fields V = V + 075 -I- s -I- , assumed to be weak. Our 
purpose is to evaluate the quark propagator S, which is 
just the inverse of the argument of Det in 

We can expand the quark propagator S with respect 
to a single instanton field Ai: 

S = So + J2{S, - So) + ^(5, - So)S-\S, -So) + . .(12) 

i i^j 

Now we are going to make expansion over V and express 
Si via Si- Notice that the field V transforms nonho- 
mogeneously. If we could evaluate Si exactly, the final 
result would be gauge-covariant. However, due to the 
zero-mode approximation expansion over V would kill 
the gauge-covariance of the propagator Si. In order to 
preserve it, we introduce auxiliary field 



= L,ip + V)L,~p, 



(13) 



and the gauge connection Li, which is defined as a path- 
ordered exponent 



(14) 



Zj(x, Zi) = 74^1(2;, Zj)74 



where denotes an instanton position. The field 
V^{x, Zi) under flavour rotation 

ip{x) — * U{x)ip{x) 

transforms according to 

Vl{x,Zi) ^ U-\zi)Vl{x,z,)U-\z,) 

and thus expansion over it does not violate the gauge 
covariance. The propagators Si and 5*0 in terms of the 
field Vi have a form: 



Si — LiSiLi, Si 



1 



p + Ai + VI + im 
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(15) 



P + V/ + im 

Expanding S- over y/ and re-summing it, we get 

si^s.,{i+Y,i-v:s,r)= (16) 

n 



where 

b,^{'^>o,\piS'o,~So)mo^), (17) 

Q - b, = -{^o^\pS'oiP\'^'Q^) = 

($o.|(«m + V:/)|$o^) - {<^>o^\i^m + V;)SUim + V:)\^o^) 

Rearrangement of the Eq. for the total propagator 
yields 



S = So + SoJ2 L-'m^o) ( ^ + ^C^^ + • • • ) (*o,|pL7'5o 



-V -T 



{^oj\pL~^So 



(18) 



where 



(<i>odp(L-i^oZ-i)p|<i>o,) - {^o^\pSoL-^L,p\^oo), 
= (1 - %)($o,:|p5oir'-^j-pl*o,), 

{l-6.,)V.,. (19) 



It is natural to introduce 



1 



|0o) --Lp|$o) 



(20) 



which has the same chiral properties as zero-mode func- 
tion |<i>o). Then 



1 



y + T 



S - So^ -'5'o^p|(/)oi)((/!)oi| 
with 

V + T = fSop. 
From pT|) we get 

Tr(^-S'o) = - ^{(^)a,]\pSlp\(j)Q,^){(|)o,, 

id 

Introducing now the operator 



|</'oj)(0oj|p^o(21) 



(22) 



/ 1 



\pSop 



(23) 



fa 



where indices i,j refer to the number of instantons and 
/, g are flavour indices, it is easy to show that 



In (Detio„) = iTr / dm'{S{m') - So{m')) = (24) 



dB(m')ii,~, ,,,1 ^ , ^, s 
■ ^ ^'■'-(S(m'))-/ =TrlnB(m) 



dm 



■1,3 

Thus we have 



Detiow [v, a,s,p,m] = det B{m) 



(25) 



where i? is the extension of Lee-Bardeen's matrix B [26[ 
for the case of nonzero external sources {v,a,s,p) and 
current quark mass m. Certainly, in the limit of small m 
and zero sources B coincides with B. 

Independent averaging of the Detiow over instanton 
collective coordinates ^, justified in the dilute gas ap- 
proximation and performed with the help of fermioniza- 
tion [l^, [lO, HH gives the desired light quarks partition 
function Zn[v, a, s,p,m\: 



Zfq = (Detiow b, a, s,p, m]) = (deti?) 
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where 



(27) 



arc the t'Hooft-likc non-local interaction term with Nf pairs of quark legs, 



(28) 



and gauge links are defined in (fT4|) . 

The form-factor of the nonlocal interaction is com- 
pletely defined by the quark zero-mode. The fermion 
fields i/j\4' s-rc interpreted as constituent quarks. For 
the exponentiation we use Stirling-like formula 



W. 



N± 



dX± exp{N± In 



A± 



N± + X±W±), 



(29) 



where X± play a role of the dynamical coupling constant, 
defined by saddle-point condition at the integral ([29|) . 
Notice that the validity of the formula (|29p is controlled 
by the large number of instantons N± . 

The partition function in the instanton vacuum after 
integration over collective coordinates and exponentia- 
tion may be reduced to the form (for the case Nf = 2): 



dX+dX_DijjDipe- 
K 



(30) 



5 = iV± In iV± + V^t(i(9 + f + im)i^ + X±Y^il) 

X± 



y± = dot J± -t- /?2 det J± 



1 



1 



8iVc 2iVc - 1 BiVc 



= o[ — 

Nr. 



(32) 
(33) 
(34) 



T± / t' 1 ± 75 , , 1/^7,* T 



SmCC cr^^ = ^''A'J'pA^pA 

L = exp (^i I V^iOd^^^ , L = joL^ 



70 



, (35) 

(36) 

(37) 
(38) 

(39) 



where the determinant is taken over the (implicit) flavour 
indices, and K is some inessential constant introduced to 
make the argument of logarithm dimensionless. From 
([33|) one can clearly see that the contribution of the ten- 
sor terms is just a 1/iVc-correction. For the sake of sim- 
plicity wc will postpone consideration of the tensor terms 
contribution until the Section rVIII and concentrate on the 
first term in ((32|) . 

Next step is the bosonization of the interaction term 
W± , ip] . Here we are considering only the Nf = 2 case, 
for which bosonization is an exact procedure, and take 
N^ = N_, as discussed before. So, we get the partition 
function 



ZN[v,a,s,p,m] 

= J dXD'i>DijDi;exp N\n^-N + J dx ^-2$^ + ij'' {p + V + im + i ^'^'^^^ ^ LF<^> ■ TFL)i' 
or, integrating over fermions, 

where, 

K f o I p + V + im + i 

-S[X,<P,v,a,s,p,m] ^N\n - 2 / (i.T$2 + Trln 

A J \ p + v + im 



2g 



-LF^ ■ TFL 



(40) 



(41) 



where 



(Af^-l)2JVe 
2No-l 



is a color factor. Here F = {1, 75, if, if75} correspond to the decomposition $ 



6 



{$0, ^} = W, ^1 $}i and = + = cr^ + r?^ + 
5"^ + (jr. Notice that in (PU)) in contrast to the NJL model 
the couphng constant A is not a parameter of the action 
but it is defined by saddle-point condition in the integral 

Note that the partition function (|40p is invariant un- 
der local flavour rotations due to the gauge links L in 
the interaction term V±j[ij}'^ , -ip]- However, instead of the 
explicit violation of the gauge symmetry due to the zero- 
mode approximation (|6l7p . we have unphysical depen- 
dence of the effective action on the choice of the path in 
the gauge link L. In our evaluations we used the simplest 
straight-line path, though there is no physical reasons 
why the other choices should be excluded. In Sec. 1x1 we 
demonstrate explicitly that for the quantities evaluated 
in this article the path dependence drops out. 



III. DYNAMICAL QUARK MASS 

One of the main advantages of the instanton vacuum 
model is the natural description of the SxSB, which is sig- 
nalled by non-zero vacuum quark condensate {qq). The 
quark-quark interaction term (^7)) leads to the strong at- 
traction in the channels with vacuum (and pion) quan- 
tum numbers. As a consequence, there appear the 
nonzero vacuum expectation a of scalar-isoscalar com- 



ponent of meson fields $ and related with it {qq). For 
evaluation of the partition function Z[m], it is very conve- 
nient to use the formalism of the effective action [ssl . [s^l 
re//[7Ti, A, $], defined as: 



(42) 



J dAcxp(-re//[m,A, $]) 



where for the sake of simplicity we dropped all the exter- 
nal currents which are not essential in this section, and 
the field $ is the solution of the vacuum equation 



9re//[m,A,^] 



0. 



(43) 



Notice that the solution depends on A, i.e. ^ = $(A). 
Up to the Section fV Al it will be assumed that the only 
nonzero vacuum field is a condensate $ = tr, which 
is independent of coordinates, so the effective action 
Teff[m,X,^] may be replaced with effective potential 
Veff[m,X,a]. 

In the leading order, the effective action just coincides 
with the action ([4T]). Shifting ^ ^ a + ^' and inte- 
grating over the fiuctuations, we get for the meson loop 
correction 



J 



r-//[m,A,a] 



^Trln ( 45,, , - ^rTr 



M{p) 



M{p) 



p + i^i{p) p + ifi{p) 



(44) 



r 



where fi{p) ~ 'm+M{p) and we introduced the dynamical 
quark mass M{p) = MF'^{p); M 



2s 



-a. 



It is convenient to introduce notations for the leading 
order meson propagators 

1 



V^{q) 



x2 "^2 



(45) 



where 

V^iq)=TriQ{p)T,Q{p + q)T,) 
and 

iM{p) 



Q{p) = S{p)iM{p) 



p + ifi{p) ' 



(47) 



In the leading order (LO) over Nc this equation simplifies 
to 

1 



4ct" 



V 



TriQip)). 



(50) 



The equations (|48l50p completely define cr, a^o as 
hmctions of A. Using ((50)) . we may cast (|45|) into the 
form 



(46) Il-\q)^n-\q)^ 



1 



{TrQ{p) + V^iq)) 



(51) 



With these notations vacuum equation (|43|) turns into fl~^(0) — 8 m A^c f 



This replacement guarantees self-consistency of the 1- 
mcson-loop approximation and complies with Goldstone 
theorem: 

df^p 



da 



Aa^-^Tr{Q{p)) 



(48) 



M{p) 

(27r)>2+^2(p) 



(52) 



2m{qq){m = 0)\lo 



1 f d^q 

W) 

where 

V;{q)^Tr{Q\p)T,Q{p + q)T. 



^^(.)n.(,) = o, 



(49) 



Notice that the pion propagator has a correct chiral pole 
if and only if the field a satisfies the LO equation (|50p . 
Similar property is observed in NJL model (see, e.g., [15[ 
and references therein). 
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There is an important difference between the instanton 
vacuum and traditional NJL-type models - the coupling 
A is not an external parameter of the model, but is de- 
fined from the saddle-point equation. We have to inte- 
grate over the coupling A in Eq. (|42)) to obtain partition 
function Z^. The saddle-point approximation for the re- 
sult becomes exact in the large- limit (see discussion of 
(P^l above). The saddle-point equation for A has a form 



A 



dX 



(53) 



E(^2('?)-^'(9))n.(9) = o 



Notice that V2-term in ((53|) requires special attention. 
Formally it is next to leading order correction, while nu- 
merically it is strongly enhanced (about a factor of 30 
compared to the other l/A^c-corrections), which indicates 
the failure of the large- iV^ expansion in ([551) . 

If we solve the equation (j53|l . expanding it in powers 
of 1/Nc with the set of parameters Q, we'll get 

(54) 



Mo = 0.567 - 2.362 m 
Ml = -^(-0.687 - 0.808 m 



4. 197 m In m). 



Here and in the following M and m are given in GeV. 

We can sec that the meson loop correction Mi is of 
comparable size with the LO term Mq, so we can try 
to solve the equations ^8\53\) numerically in chiral limit 
and then evaluate the chiral corrections to it. Such pro- 
cedure gives 



A/(m) = 0.36 



Tfi 

2.36m- — (0.. 



4.1971nm) (55) 



The accuracy of the solutions (|54I55[) is 0{r 



1 iv? 



FigQ] 



represents the M(m)-dependence obtained from Eqs. 
(|48|53p . For the sake of comparison we also plotted the 
lattice data from [s^ . From the first point of view our 
results are w 30% higher than the lattice data. However, 
they are given in different gauges and on different scales. 
Since M{p) is essentially nonpcrturbativc object, it is 
not very easy to rescale the data and make comparison. 
Rough estimates in perturbative QCD show that the dis- 
crepancy may be attributed to the scales difference. We 
may conclude that we have a qualitative correspondence 
between our model result for A:/ (m)-dependence and un- 
quenched lattice data [35[, as it was expected. 



IV. QUARK CONDENSATE 

The presence of the quark condensate (qq) is one of 
the most important properties of the QCD vacuum. Its 
value characterizes the SxSB. In the chosen framework 
we can extract it directly from the effective action taking 
derivative over the current quark mass [36| 



1 d^eff 

2 dm 



.-iTr 
2 



p + in{p) 



I Bmvmciii, imijiieiidwit 



M(m) 



FIG. 1: m-dependence of the dynamical quark mass M. The 
solid curve - the exact numerical solution (|55|) of vacuum Eqs. 
(|48l I53p . The dashed curve - the solution (|54|l . obtained by 
the iterations (1/A^c-expansion) with the same accuracy. Data 
points are from [s^ . Notice that the scale of the lattice data 
is 1.64 Gel/, not p"^ « 0.6 GeV. 



P- 



d'^q 



i2n) 



^l/««(<z)H,(g), 



(56) 



where the first term is a contribution of the tree-level 
action and the corresponding vertex part V^'^{q) is 



= Tr 



MF^p) MF^{p + q) 

'p + q + i^i{p + q) 



Evaluation of (|56)) gives 

-{qq){m) = ((0.00497 - 0.0343 m) TV^ 
+ (0.00168 - 0.0494 m 

-0.0580 m Inm)) [GeV^] 



(57) 



The (^q)(m)-depcndence is depicted on the Fig. [2] We 
can see again that due to the chiral logarithm the m- 
dependence is not linear and meson loops change the 
m— dependence of the quark condensate drastically. 

Here and below we define as leading order (LO) the re- 
sult calculated with the l/A^c-corrections dropped every- 
where except in the mass Af(m), which is taken from (j55p 
without the last ©(m/A^c, "^/^c lnm.)-terms. As ncxt- 
to-lcading order (NLO) we define the contribution of the 
last two terms in (|55p plus " direct" contribution of meson 
loops. 

For the sake of comparison, we also plotted the value 
{qq)o which one would get using LO formulae with the 
mass Mo{m) taken from (jM]) . Recall that the value 
{qq(rn = 0)) = (255 Mey)^, as weh as F^{m = 0) = 
88 MeV, was used as the input in order to fix the parame- 
ters (p, i?) in (21). In LO we have {qq){m ~ 0.1) /{qq){m = 
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FIG. 2: m-dependence of the quark condensate —{qq). The 
long-dashed curve is the LO result —{qq} lo, the short-dashed 
curve is the NLO contribution —{qq)NLO, the solid curve is 
the total contribution —{qq) lo+nlo- The dot-dashed line 
represents the leading-order in 1/A^c-expansion result, evalu- 
ated with the mass Mo from (|54[) (see text). 



0)|lo = 0.31; while with NLO corrections we have 
{qq){m = 0.1)/ {qq){m = 0) = 0.89, a noticeably different 
result showing the size of the chiral logarithmic term in 
the NLO-corrcctions. 



V. QUARKS IN EXTERNAL AXIAL- VECTOR 
FIELD AND PION PROPERTIES 



The formalism of the effective action used in one of 
the previous sections may be successfully applied in the 
presence of the external axial- vector isovector field = 
The general partition function (|1D|) is reduced in 
this case to the form 

Zjv[m,a^] = J dXexp{-Teff[m,\,u,afj]) (58) 

The external field can generate nonzero vacuum aver- 
age {(f)) = u and shift the value of the vacuum filed (cr) 
and saddle-point value A. 

In this paper we restrict ourselves to the case of the 
soft and weak external field 0^(5), which can be treated 
in perturbative fashion, and q ~ M^r <^ p^^. For 
the purpose of this paper it is sufficient to keep only 
0(5^, a^9pU, (9pM(9pM)-terms in ([55]). 

On general grounds, one can guess that the vacuum 
expectation value m a^, whereas shifts of the vacuum 
field (ct) and saddle-point value A are proportional to the 
second power of a^. Using the saddle-point equation 



dTeff[m,X,u,d^ 







(59) 



and the vacuum equations 

dTeff[m,\,u,ai,] ^ ^ are//[m. A, a^] ^ ^ 
da ' du ' 

we may easily get that the shifts of (cr) , A contribute 
only to C'(a^)-terms and thus may be safely omitted in 
this paper. 

In the leading order the effective action simply coin- 
cides with the classical 5* defined in (|4T|) . Using vacuum 
equations (|60p . one may show that (ct)^ -I- ((/))^ = const. 
This inspires us to introduce a unitary matrix U with the 
properties 



U = uo + iTu, U'^U^UU'^ 
{a) = (Tito, ~ fit. 



(61) 
(62) 



where a is the value found in Section IIIII In this rep- 
resentation the vacuum meson field is represented as 

In the next to leading order one has to take into ac- 
count the fiuctuations of the field <I> all + and 
integrate over <&'. 

Meson loop contribution to T^f / has a form 

(5^5[m, A, a, u, a^, $'] 



1 



rr/;[m,A,i7,a^] = -Trln 



*'=0 



= r™// [m. A, M = 0, = 0] + Ar;7/ [m. A, u, a^] (63) 

It is convenient to rewrite the second derivative 
5'^S/5^'^ as 



J 



(64) 



where Ai^ij is of the first order in external field a and 
induced field u, and ^2,ij is of the second order in them. 
Then 



1 



1, 



^prnes ^ Trn,^^^.^^,,:, " -Tr[n,<5„- Ai,,,] 



(65) 



Expanding (j65p and collecting the terms a^a^, a^j^d^Ui 
and d^Uidfj^Uj, after simple but very tedious evaluations 
it is possible to show that in agreement with chiral sym- 
metry expectations, the structure of the effective action 
is 

Teff = ao{dfj, + di_,uf + m aidf_,u{af^ + df_,u) + (66) 



+ma2U ^raa-p 

where the constants Fij differ only beyond chiral limit: 
P'L - Fun = 2 {Fl - i^i) = -a, m. (67) 



From (|66l67p one can get that the two-point axial- 
isovector currents correlator has a form: 



2 I e It^Qi^ 



(68) 



= 5,,F^ S^. - 



q- 



Mi 
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We can see that M^^ has a meaning of pion mass and 
- pion decay constant. Numerically it is much eas- 
ier to calculate F'^ as a constant in front of (J^^-tcrm in 
taking u{x) = 0, a^{x) = const (this corresponds to 
a^(g K, 0)). Also, it is possible to show that the result of 
such evaluation is independent of the path choice in the 
transporter L (see Section |X] for more details). In a simi- 
lar way, we can put a = Q,u — u{q) and evaluate in NLO 
the quantities F^^ and pion mass M^. Both quantities 
F^ and naturally have the chiral log terms due to the 
pion loops contributions. The coefhcients in these chiral 
log terms are controlled by the low-energy theorems [l^l 
and are reproduced analytically in Section [IXI 

We can see that the specific structure of the Tef / ([66]) 
provides a check of the numerical calculations. Moreover, 
the chiral log theorems provide another check of the nu- 
merical calculations. 



— mlnm + 0(m2) •10"3[Gey^]= (69) 

(7.67- 11.35 m- 44.25 m In m) • 10"^ [GeF^] 

TO is given in GeV and the constant in front of 0{mP) 
is given in GeV^, the constant in front of 0{m) is given 
in GeV. The i^7r("i)-dependence is shown in the FiglS] 
For the sake of comparison, we also plotted the value 
-Ftt.o which one would get using LO formulae with the 
mass AIo{m) taken from ((54|) . Recall that the value 
F„{m = 0) = 88MeV, as well as {qq{m = 0)) = 
{255 MeV)^, was used as the input in order to fix the 
parameters {p,R) in The comparison between the 
solid curve and the long-dashed one shows that the effect 
of the NLO-corrections grows with to and is about 40% 
at m = 0.1 GeV. 



A. Pion decay constant from a'^-term 

The basic diagrams which contribute to this quantity 
in the leading order and in the next-to-leading order are 
shown schematically in the Fig. [3] and Fig. 3] respec- 
tively. Notice that in integration over </> the saddle-point 
is shifted to {(j)) ^ a^, where the "proportionality" sign 
implies some nonlocal linear operator. All the vertices on 
these plots should be understood as a sum of the local 
and nonlocal parts, as it is explained in the Appendix \X\ 
and in the Fig. [TUl 




FIG. 3: The basic diagrams which contribute to the term in 
Fe//. The wavy line corresponds to the external field a^{x), 
the dashed line corresponds to the intermediate meson, the 
bulbs correspond to all the possible (local and nonlocal) cou- 
plings of the field a to the constituent quarks. 




FIG. 4: The meson loop corrections to the term. The 
notations are the same as in Fig. [3] 

Finally we have 



F jm), GeV 



0.2 _ F,,, 




FIG. 5: m-dependence of the pion decay constant Ftt. The 
long-dashed curve is the LO contribution, the short-dashed 
curve is the NLO contribution, the solid curve is the total 
LO-I-NLO contribution. The dot-dashed line represents the 
leading-order in 1/A^c-expansion result, evaluated with the 
mass Afo from (|54[) (see text). 



B. M,r from pion propagator 

Effective action Fg//, Eq. (|66|) at a = has a mean- 
ing of inverse 7r-meson propagator at small external mo- 
mentum q Afjr with account of meson loops. For our 
purpose it is sufficient to have only the 0(9*^) and 0{q^) 
terms. 

The LO propagators of the mesons were defined in (|45|1 . 
The inverse pion propagator with account of meson loops 
is given by 



10 



(T 



1 



ndm,{k + q)Vf{k,q)Vf{k,q) 



a- J (27r)4 '^^ ^ '-^^ ^6 J (27r)4 

where the function Q{p) was defined in (|47|) . and the vertices are 

V^{k, q) = 2Tr (Q(p)i75T'^Q(p + g)r.Q(p + g + k)TjQ{p + q)iToT') 
Tr (Q(p)z75T'^Q(p + 9)r.Q(p + <Z + fc)*75r''0(p + fc)r,) 
F/(fc, q) = Tr (gb)»75r'^Q(p + q)T,Q(p + q + k)V,) 



(70) 

(71) 
(72) 



FIG. 6; Diagrams corresponding to the fast two meson ioop terms in (|70|) 



The value of the pion mass Mjr is defined as a pole 
position in the propagator ([70)) . Now we would like to 
demonstrate that ~ m, i.e. that the loop corrections 
do not violate the Goldstone theorem, i.e. 11^^(0) to. 

In the leading order over the only contribution to 
Ii'^^[q) comes from the first term in (|70p . and using equa- 
tion ([5D|) may be reduced to 



n-i(o) = TO 



d^p MF^{p) 

(27r)4 p2+^2(p) 



(73) 



In the next-to-leading order contribution of the first term 
in (iTOl) is 



1 



-n,(g)V3(g) -Hto- 



(27r)4 p2+^2(p) 

d'^p MF^{p){p^ - APf{p) + to2) 



(27r)4 '^^^ ^^^^ Mo \J (27r)4 

and explicit expression for V^{q) defined in ((49l) is 



(p2-f;,2(p))2 



(2^ 



(74) 



APF^p)F^p + g) (±2 (p^ + p ■ fc) ^(p) - p^fijp + q)+ m'(p)m(p + g)) 

(p2+^2(p))2 ((p + g)2^^(p + ^)2) 



where here and below 

1 , i = (T, 77 
—3, i ~ (7,(j) 



(75) 



77, (/) and — for a, a 

The vertex part V^(k,q = 0) in the second term of 
([70)) may be rewritten explicitly as 



and the sign in the numerator of the first term is -f for 



11 



(76) 



(2^ 



2c, 



Ad^F^{p)F^{p + k){fi{p)n{p + k)±{p'^ +p-k)) W^F^ip) M'^F'^{p + k) 



I 



The third term in ([70]) may be reduced to the form 



d^k 



(n,(fc)fir(fc) + n^(fc)nff(fc))y2(fc) (77) 



(27r)4 '"^'^ 0V 
where F(fc) equals 

d'^p M^F'^ip) M F^ip+k)^i{p + k) 



(27r)4 p2 + ^2(p) ^ fc)2 + ^2(p _^ fc) 



(78) 



To ehminate the terms with double propa gato rs in the 
third term, we follow the strategy used in [37| for pure 
NJL. Using explicit expressions for propagators, we can 
notice that 



l6Nc f d^p MF^{p)fi{p) MF^{p + q)fi{p + q) 2V{q) 



where 



(27r)4 p'^ + fi'^{p) {p + qf + iJ?{p + q) 

vn=8/V [ "^^P MF^{p)fi{p)MF^{p + q)fi{p + q) 
'J (2^)4 p2+^2(p) (p + qy+^^p + q) 

Obviously 

V{q) = V{q)+m6V{q) 



(80) 



SViq) 



Hence 



-8N,. 



d^p MF^ip) MF^{p + q)n{p + q) 
(27r)4 p2 + ^2(p) (p ^ q)2 _^ ^^2(p ^ 



fl,(g)fl^-(q) = 
n^(g)n5(g) = -a 



n^(g)-n^-(g) 

2V^(g) 



2V{q) 



(82) 
(83) 



(79) 



and (|77p turns into 



(2^: 



4 (n,(fc) + n,(fc) - n^(fc) - n^-(fc))F(fc) x 



(84) 



We can notice that coefficients (signs) in front of each of 

2 



the terms in (|84p may be rewritten as ^^ii . Finally after 



(81) several simplifications we arrive to 
I 



8Nc 



d'^k 



U,{k) 



d^p 
(2^ 



2c, 



MF^{p) MF'^{p + k) 



{fi{p)fi{p + k)±{p'^+p- k))M^F'^{p)F^{p + k) 

ip' + fi^{p)f [ip + kf + ^lip + k)^) 

2 (ci + 1) 



(85) 



(m + Cip{p)) 



p^ + ^^^{p) {p + k)^ + ^'^{p + k) 
8Nr. Ml ( f d^p MF^{p)ip^ - APf^ip) + m^) 



(27r) 



V{k) 



il/n 



(27r)4 



{p^ + fi^p)y 



We can see that the pion propagator 117^ (0) ~ m, i.e. the 
loop corrections do not violate the Goldstone theorem. 
Notice that the function V{k) is positively defined, so 
no unexpected poles should appear due to the term oc 
f jtkU fk)^^^ 

J (27)^-^-^^'^'' Vik) ■ 



Evaluating i^^„-term in ([5^ with account of meson 
loops and taking the ratio of ([55]) to F^^, we obtain the 



pion mass 



Mi = m 3.49 



1.63 

18.25 13.5577, , 
m [ 15.5 H 1 Inm + (^(m^) 



(86) 



m(4.04 +21.587TO + 4.52mlnm + 0(7n2))[Gey2] 

in the next-to leading order. The Af^ (TO)-dependence 
of the pion mass is shown on the FiglT) For the sake 
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of comparison, we also plotted the value M-^^ which one 
would get using LO formulae with the mass MQ{m) taken 
from ([54|l . Altogether, for this observables the NLO- 
corrections turn out to be small. 




m, GeV 



0.025 



0.05 



0.075 



0.1 



FIG. 7: m-dependence of the pion mass Mt^- The long-dashed 
curve is the LO contribution, the short-dashed curve is the 
NLO contribution, the solid curve is the total LO-I-NLO con- 
tribution. The dot-dashed line represents the leading-order 
in l/A*'c-expansion result, evaluated with the mass Mo from 
((54)) (see text). 



VI. FINITE WIDTH CORRECTIONS 

In the previous sections it was assumed that the in- 
stanton size distribution has a zero width, i.e. 



d{p) = 5{p - p) 



(87) 



and all instantons have the same size p. As it was shown 
in the early works of Diakonov, Petrov^, [l3|, this ap- 
proximation is justified by the small parameter 1/Nc, i.e. 



{P? 



For numerical evaluations we take the value 



5p' = ip') {P? 



0.5599 GeF- 



(88) 



(89) 



which follows from the two-loop size distribution of Q 
Since we are interested in all \/Nc corrections, we must 
take the finite width into account. To do this, we must 
return to the formula ((28)) . Additional integration over p 
doesn't change the exponentiation procedure, and we get 
the standard 27V^-interaction term in the effective action 
S. However, for bosonization we should slightly modify 



the standard procedure. Here we consider only the case 
Nf = 2 we are mainly interested in. For this case 

J D^{z,p) exp^-iy d'^zdpd{p)<^'^{z,p)+ 
d'^zdpd{p)^{z, p)J{z, p) 



where $(0,/o) are pure mathematical objects (matrices) 
Thus the effective action turns into 



5=^lnA 



2 j d'^zdpdip)^^iz,p) - (91) 

Tr In (^p + im + ic A°-^ J d^zdp d{p)K{z, p)^{z, p)^ 

where c is some inessential constant and Kx,y{z,p) ^ 
4>{x — z, p)(t>{z — y, p). The LO vacuum equations which 
follow from ([5T|) are 

HO.5 



ic \°-^ J d^zdp d{p)<^>{p)K{z, p) 



N 1 / 
— = -Tr - 

V 2 \p + im + icX^-^ ^ d^zdpd{p)K{z,p)<^{p) 



(92) 



4$o(p) - Tr 



icX°-^K{z,p) 



p + im + ic A"-^ / d-^zdpd{p)K{z, p)^o{p) 



where $o(p) the quantum numbers of a. The un- 
known variables in these equations are parameter A and 
function ^{p). In general case these equations are very 
complicated and may be solved only numerically. How- 
ever, for the case of small width we can make expansion 
over Sp^ = {p'^) — (p)^. In particular, we expand the func- 
tions ^{p) and K{p) over (p — p), get a system of equa- 
tions for the Taylor coefficients $i and A. Solving these 
equations, we can restore the function $(/o) (at least in 
the small vicinity of p) , and after that evaluate the finite 
width contributions to all the quantities. Details of these 
evaluations are given in Appendix |B] The final results 
for the corrections are 



Sp^ 3 

(0.0024 - 0.020m + Omdin^) [GeV^] 
Sp^ 3 

(-0.0024 + 0.058m - O.SSm^) [GeV^ 



Substituting the numerical value ((89)) . we get 

(0.00045 - 0.0037?n + O.OOSGm^) [GeV^] 
S{qq) = 

(-0.00045 + 0.011m, - 0.062m2) [GeV^] 



(93) 
(94) 

(95) 
(96) 



Here m is given in GeV. Thus we can see that these cor- 
rections are relatively small, w 5% for and « 2.6% for 

im)- 
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If we define the mass M{p) of the constituent quark as especially simple: 



M{p)^cX'>-^ J d^zdpd{p)K{z,p)<i>o{p), 
M{p) = MF'\p) + 5p^5Ah{p) + O {5p^) 



(97) 
(98) 



we can get the p-dependence shown in the Fig. [HI We can 
see that for p = Q the increase of the dynamical quark 
mass is 5M/M = 10%. 



M{p), GeV 



M(p)+5M(p) 

m(p) 

M(p) 




_l I I I I I I L_ 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 p. GeV 



FIG. 8; Change of the p-dependence of the constituent quark 
mass M{p) due to the finite width corrections. The dashed 
curve is the contribution of the leading order result, the dot- 
dashed curve is the contribution of the finite width correction, 
the solid curve is a total result 



VII. TENSOR TERMS CONTRIBUTIONS 

In this section we study contribution of the {1/Nc- 
suppressed) tensor terms to the axial currents correla- 
tor (a^Oiy) introduced in (|30|l . Notice that due to the 
identity ((551) f^e fields Jj^'^ depend on each other. Under 
Lorentz transformations J^^i/ transforms like (1, 0)-|-(0, 1) 
whereas chiral components J^^^ transform like (1, 0) and 
(0,1). 

In case Nf = 2 evaluation of the determinant becomes 



det J± - - I Vfj(*))2 _t 



F, = {1,zt} 



(100) 

(101) 
(102) 
(103) 



The term ±^eap^ujjiv jj^p effectively averages to zero if 
the vacuum is symmetric, 7V_|_ = For bosonization 

we use standard formula 



dye 



-y —axy 



Notice however one subtle issue which exists for the ten- 
sor fields. The field J^i, is antisymmetric, J^^ = —Jfiv, 
thus not all the components are independent and this 
fact must be taken into account in bosonization to avoid 
double counting : 



exp(a2j2,/8) = (105) 
/ D^f,^ cxp - E - a E J^..^-^,. = 

= y" i?$^.exp - 

The above-mentioned antisymmetry must be also taken 
into account in all further differentiations over ^^jy. 
From (|106|) we get for the partition function 



d\j^d\^D^Di;D<^^D^%e 



5 = -iV±lnA± +2 



i$2 

2 2,/J,Z> 



iX'^ '^LFip) (^a$,F, + -/3$,,^.a^.F, ) F{p)L- 
We can integrate out fermions and get 
Zn[v,T]^ J d\dXD^D<Pf,^e-^ 



S = -iV±lnA± +2 ( $2 
Tr log 



1$^ 

2 2,/J,Z> 



(106) 
(107) 

(108) 
(109) 



p + v + -T^^o-^^ + im- 



1 



A°-5ZF(p) a$,F, -I- -/3$,,^.a^.F, F{p)L-^ 



where we have introduced external tensor source T^i,. 
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FIG. 9: Contribution to the axial currents correlator. The 
intermediate state is the tensor meson ^^v, the diagram is 
1 / A^c-correction. 



1 - 



aia 



(2^)4 (p2+^2(p))2 



Thus we can see that this diagram is just 0{q^)- 
correction to the axial correlator. Since we are interested 
only in the LO over (evaluation of F^r), we should not 
evaluate this diagram. 



VIII. GASSER-LEUTWYLER COUPLINGS 



A. Contribution to the axial currents correlator 

In this section we consider the contribution of the ten- 
sor mesons to the axial currents correlator. This 
contribution may be represented as a Feynman diagram 
shown in the Figl^l Straightforward evaluation of the 
tensor-axial coupling is 

'^iifj.iypx^fiiyqpaxiq) x ca, (110) 
d^p 2ti\p)+pMf{p)r{p){p^ - 3/x2(p)) 



According to [l|, the low-energy constants U of the 
chiral lagrangian may be extracted from the 0(m)- 
corrcctions to physical quantities, e.g. 



Mi 



CA 



327r2F2 

,2 



m 



1 + + 



(112) 
(113) 



(2^)4 (p^ + fiHpW 

and the total contribution of the diagram in the Fig. [5] is 



(111) 



J 



where Mt^^F^^ are the pion mass and decay constants, 
= 2mB and B^F are the phenomenological param- 
eters of the chiral lagrangian. Using our results 
we can obtain (m is given in GeV) 



= 0.00284777A^c - 0.000868917 -f O ( 



(114) 
(115) 



-1.14251 iVr 1 



0.872354 I 0.874555 In m 



N, 



0.936972 



o 



15 In m I / 1 

' 1" '-^ I IV^ i i /I 

^ - 0.0738267- 1.14251 TVc- 0.999 Inm + O ( — )(116) 



-0.0793814 A^, 



232149 I 12.5977 In m 



0.468486 



o 



(*) 



-0.0793814 7Ve + 0.0187608- 1.000 In m + O ( — ) (117) 



which gives 

F = 88MeV, B ^ 2.019 Ge^, I3 = 1.84, U = 4.98 (118) 

at m = 0.0055 GeV, corresponding AU = 0.142 GeV, 
F^ = 0.0937 GeV. 

The values of F, -{qq{m = 0)) = -F^B in (fTT^ were 
taken as input when we fixed the parameters p, i? in ([4]). 
Our values of (^3,^4) should be compared with the phe- 
nomenological estimates [12, H^l as well as lattice predic- 
tions [H, [3^ given in Table |T1 



Notice that in (fTTBl - [TT^ we keep 0{Nc, N°)-teTms 
and drop NNLO terms 0{1/Nc) in agreement with our 
general framework. Without such expansion one would 
get 

[3 = 0.28 h = 4.28 (119) 

IX. CHIRAL LOGS 

Now we would like to discuss the chiral logarithms 
TO In m generated by pion loops in the effective action 



15 





XPT 


MILC 


Del Debbio 


ETM 


Our 




[1, 23,40] 


m 


et. al. [39] 


[42j 


prediction 


^3 


2.9 ± 2.4 


0.6 ± 1.2 


3.0 ±0.5 


3.62 ±0.12 


1.84 


k 


4.4 ±0.2 


3.9 ±0.5 




4.52 ±0.06 


4.98 



So evaluation of the chiral logs is straightforward but 
tedious. 

The first relation for the constituent quark mass M (Tsj 



TABLE 1: Estimates and predictions of the low-energy con- 
stants. The first column contains phenomenological esti- 
mates, the next three columns are lattice results from different 
collaborations, the last column contains our results. The first 
four columns of the table are taken from [231 . 



M(m) = Mo ( 1 - 32^5^^ l^^'^'.^o + 0(A4') j (122) 



Teff- First of all, we should notice that we are mak- 
ing not just a chiral expansion, but a double expansion 
over the 1/Nc in large- iVc limit and the chiral expansion 
over the quark mass m. Wc must take this fact into ac- 
count when checking the thcorcms-wc should keep only 
LO and the NLO corrections. Since chiral logs originate 
from pion loops, which arc C'(l/A^c)-corrections, the gen- 
eral chiral log theorem in our approximation has a form 



G{M^)=G{AU = Q) 1±7 



where 7 is some specific numerical coefficient, and we 
took into account that j^q ^ Nc- 
Chiral logs have two sources: 



1. Shift of dynamical quark mass 5 Ad due to meson 
loops (" indirect" meson loops contribution) . In this 
region contribution to chiral logs is obtained acting 
with operator ^gfy on the analytic LO expres- 
sions for the correlators. 

2. Pion small-momentum g « region in pion loops 
("direct" contribution). In this region pion propa- 
gator may be approximated as 



may be most easily obtained expanding the vacuum equa- 
tion ((53l) and collecting 1/A^c-terms. 



(123) 



\ J (04E(^2(9)-^3('7))n.('z) = o 



Using (|120p and the relations 



M2 



(2^)4 q^ + Ml^o 



Mn 



d 



Tr[Q(p)] = 16iV, 



dMo 

y/(o) -y/(o) = -8iv, 



167r2 --^'^o 

(2^)4(p2 + M2(p))2 
{27r)Hp^ + f^^pW 



we can check that chiral log in Mi satisfies p22p . 

In completely the same way, expanding Eqs. (j56p and 
using (|122p . we can obtain for the chiral log in quark 
condensate 



LO 



vertices Auj and A2^ij 
dent of g), so 



(120) 

as constants (indepen- 



{qq{m)) = {qq)a 1 



d^a MwnlnM^ 

n(g)A(<z) « A,{q = 0)^hL2^hLO +(.121) 



(277 



167r2 



J 



Similar evaluation of the chiral logs in F^^ and Mjr 
yields 



F^im) = F'iO)ll- 



InM 



ir,LO 



0(m4) 



MlLoi^) = M 



Tr,LO 



(125) 
(126) 



in correspondence with low-energy theorems |24l | and 
large- TVc expansion. From p25ll26p we can immediately 
see that the chiral logs in the low-energy couplings k are 



correct [T|. 
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X. PATH INDEPENDENCE 

In our evaluations for the sake of convenience we chose 
the straight-hne path in the transporter Now we 

would like to demonstrate that the results of our evalua- 
tions are path independent, though on the intermediate 
steps we have explicit path dependence as a consequence 
of the zero-mode approximation ([7|). 
Let us consider first the difference of the path inte- 
grals over two contours: 6 (/ v^d£,^) = f Vf^d^^^. Due to 
Stock's theorem, this integral is reducible to the surface 
integral 



dSfiiiV^i, ix) 



(127) 



where for the monochromatic field Wp(^) = w^i((7)e''?'^ we 
have 



(128) 
(129) 



Since we are interested only in ©(g") -terms, the field 
strength is zero and thus integral over any closed con- 
tour is zero 



Vfj^d^f, ~ v^i^iq) / dS^t, = Vf,i,{q)S^,y = 



(130) 



For evaluation of Ft, from the (j^-'jjf '^)-corrclator we 
can take ~ const, which yields the path-independent 
result. 

Notice that the path independence exists only for small 
q, whereas for arbitrary q (e.g. for the "dynamical" mag- 
netic susceptibility considered in [i^) the path depen- 
dence exists and there is no other argument except sim- 
plicity why one should choose the straight-line path in 
the transporter (fT4)) . 



XI. DISCUSSION 

The aim of our work was the study of the pion 
physics beyond the chiral limit in the framework of 
the instanton vacuum model. We found the gener- 
ating functional of the hadronic correlators with ac- 
count of 0{1/Nc, m, m/Nc, m lnm/A*"c)-corrections and 
exploited it for evaluation of the corrections to different 
physical observables. The corrections considered in this 
paper include meson loops, finite width of instanton size 
distribution and quark-quark tensor interactions term. 
In contrast to the expectations, we found that numeri- 
cally the 1/A'^c-corrections to dynamical quark mass are 
large and mostly come from meson loops. As a conse- 
quence, we have large 1/iVc-corrections to all the other 
quantities. To provide the values of Fjrim = 0), {qq(m = 
0)) in agreement with xPT, we offer a new set of param- 
eters p, i? m Remarkably, this set of parameters is still 
in agreement with current phenomenological and lattice 
estimates (|2I3|) . 



The main result of this paper is the evaluation of 
the FT,{m) and Af7r('7i)-dependence with account of 
0(l/iVc, m, jTi/A^c, "^/-^c lnTO)-corrections. From com- 
parison with xPT we extract the values of the low en- 
ergy constants Our results for the values of and 
^4 are in a satisfactory agreement with phenomenological 
as well as lattice estimates (See Table H]). This means 
that the instanton vacuum is applicable for understand- 
ing of the low-energy physics, at least on the qualitative 
level. Evaluation of the other LEC's is in progress. 
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APPENDIX A: NONLOCAL VERTICES 

One of the important features of the chiral quark model 
is the nonlocal interaction of the quarks with mesons in 
On the one hand, due to the formfactors all the loop 
integrals in the model are convergent with an effective 
cut-off at the scale p ~ p~^. On the other hand, such 
regularization explicitly violates gauge invariance. The 
modification of the quark-meson coupling to 



L(x, z)F{x - z)$,(z)r,L(z, y)F{z - y) 



(Al) 



where L- factors are defined in (|M1) , formally restores the 
gauge covariance. However, it introduces unphysical path 
dependence of all correlators with vector and axial cur- 
rents. In some of the simplest cases (e.g., constant field 
a^) it is a trivial matter to show that the path depen- 
dence drops out, but it doesn't in general case. In our 
evaluations we choose conventional straight line path. 
Let us consider for definiteness the fields of the plane 
wave form 



a^(x) = a^{q)e"i'' 



(A2) 
(A3) 



The integral along the straight-line path for such fields 
equals 



d£.,v,{0 = (.^ - ;^)^«^(g)e'^- ''''^7^^ I (A4) 

iq- [x — z) 



and the same for axial field a. In general the vertex (|Aip 
contains an infinite series of the terms v, a. However 
in the next sections we'll need the vertices only up to 
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0[f ^, a^J-tcrms and restrict ourselves with this accuracy 
in this section. Then the L-factors may be rewritten as 



V^(9l)K(g2)«(2: - y)i,i{z - y)^ x 



L{z,y) ^\ + V^{q)i{z-y) 

igi -{z-y) iq2 ■ (z - y) 

L{x,z) = 7oL(z,a:)^7o = 

1 + - z)^- 7 ^ + 

iq ■ [X — z) 

V^(9i)^i'(92)«(2: - - z)^ x 



(A5) 



giiji-(a::-z) _ 2 gig2-(a;-2) _ ^ 

iqi ■ (x ~ z) iq2 ■ {x — z) 



(A6) 



where V^{q) = v^{q) + a^(g)75, V'^(9) = ^^(9) " «a'(9)75- 
The interaction vertices may be compactly rewritten in 
terms of the functions 

^iq-(x-z) _ ^ 

Gfj,{x ~ z;q) -.^ i{x ~ z)fj^———j- —F{x-z) (A7) 



iq ■ {x — z) 



H^^{x - z; gi, 92) := i{x - z)^i{x - z)^ x 

gi(ji-(x-z) _ 2 gig2-(2;-z) 



(A8) 



iqi ■ [x — z) iq2 ■ (x ~ z) 



■F{x ~ z) 



In p-space Fourier expansions of these functions are 



G,,{p-q) / d''xG^{x;q)e'^-- 



00 ^ 

00 ^ 

E (n+DKfc + l)!^^'''^^^-''"''^^-"''^^^ 



(A9) 



(AlO) 



In terms of these functions interaction terms may be 
rewritten as {^i{x) = $i(fc)e*'^^, integration over the 
dummy momenta fc, gi, (72 is implied) 



J 



L(x, z)^^(x - z)$,(z)r,L(z, y)^^(z - y) = 

[G,,(p, g)r,^^(p + k + q)+ F{p)T,,G^{p + k + q, ~q)]v^iq)Mk) + 
Gi_i,{p, qi)TiG,y{p + k + qi + q2, -(?2) + 

^ {H^^{p, qi,q2)T,F{p+ k + qi + 92) + F {p)V ,H ^,^{p + k + qi + q2] -qi, -92)) 
[-G^(p, g)r,F(p + fc + g) + F(p)r,G^(]3 + + g, -~q)]a^{q)j5^,{k) + 
- Gf^ip, qi)TiG^(j) + k + qi + q2, -92)+ 

(ifp^(p, 91,92 )riF(p + fc + gi +92) + F{p)TiH^^{p + k + qi + 92; -91, -92)) 



(All) 



W/x(9i)''^i/(92)*j(fc) 



(A12) 



a^t{ql)a^{q2)^^ik) 



+0{av) — mixing terms we are not interested in now (A13) 

Due to chiral symmetry breaking field a has nonzero vacuum expectation value (VEV). This generates momentum 
dependent constituent quark mass fi{p) and also generates new nonlocal vertices we should take into account in 
expansions over the fields $ and currents v, a: 

K - avac{[Gf,{p, q)F{p + q)+ F{p)G^{p + q, -q)]v^{q) + (A14) 
G^,{p,qi)G^{p + qi + 92, -92) + 



^ {H^i^ip, qi,q2)F{p + qi +92) + F{p)H ^^.{p + qi + 92; -91, -92)) 
[-G^(p, q)F{p + 9) + F{p)G^{p + 9, -9)]a,.(9)75 + 
- Gp(p, qi)Gy{p + 91 + 92, -92) + 

^ {Hf^i^ip, 9i,92)F(p + 9i +92) + F{p)H^j,^{p + qi + 92; -91, -92)) 



Vf,{qi)v^{q2) + 



(A15) 



aA'('?i)«"^('?2)} 



The vertices which are generated arc shown on the Fig llOl 
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FIG. 10: The basic vertices (Up to 0(^;^, a^)-terms). The 
wavy hnes correspond to the currents Vfj,,a^, the dashed hne 
corresponds to the meson, the bulbs correspond to all the 
possible (local and nonlocal) couplings of the current V to 
the constituent quarks. 



J 



$o(/o) := 00 + (f>i{p - p) + 02 



Kip, p) Koip) + Ki{p){p ^P)+ Mp) 
dpd{p) = 1 

pdpd{p) = p 

ip~pfdpdip)=S' 



(P - Pf 



N 



dpd{p)K{p, p)Mp) = <t>Mp) + y (20i/^i(p) + 4>oK2{p) + 4>2Mp)) 



-Tr. 



P,{/,c} 



^ 01 = -Trp,{f^c} 



^2 = 47^^P,{/x} 







- im + ic X^-^{KQ{p)(f)o 


+ f (20iXi(p) + 


00 ^'2 (p) + 


02ifo(rt)) 






icX^'^Kiip) 








-1171 + icX°-^{Ko{p)(l)o 


+ f (201^1 + 


00-^^2 b) + 


02i^o(p))) 






icX°'^K2{p) 








- im + ic A°-5(A'o(p)0o 


+ f (20iXi(p) + 


00 -^'2 (P) + 


02ifo(p))) 



APPENDIX B: FINITE WIDTH CORRECTIONS 



In this section we give solution of the vacuum equa- 
tions (|92p in the small-width approximation. In this ap- 
proximation we can make a systematic expansion to find 
the unknown function $(/o) and parameter A. We do not 
specify the shape of d{p), since all the finite width correc- 
tions will depend only on the width Sp^ defined in 
Expansion over Sp'^ is straightforward : 



(Bl) 
(B2) 
(B3) 
(B4) 
(B5) 
(B6) 
(B7) 

(B8) 
(B9) 
(BIO) 



Equations (jBGIBSIIBlOp form a system of four equations with four independent variables 00 — 02, ^- By definition, 



Note that none of these variables is suppressed over 6. We can make the next step and expand over the variable S: 



A:=5:A, 



5" 



;lAl 
'2A0 



From the normalization condition we get 

2 N 

20QQ = — In agreement with d = case 

001 =0 

000002 + 0?o + 000 020 = 



(Bll) 
(B12) 
(B13) 



(B14) 

(B15) 
(B16) 
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Equations (|B8IIB10|) give 

^ 'Poo = 7^^p,{/,c} ( z '. — ~~ — I -^'^ agreement with 6 = case (B17) 



4 ' ' \p + im + icloKo(j))(/)t 



00 , 



2 



^ Zi = since expression inside the brackets is not zero =^ Ai = 

-^A. ( ichKojp) \ , , . 

^ 002 = tT»^p,{/,c} -—. , . , ~ , ,^ + (B19) 



^p + im + icloKoip)4'o 
■ loKo{p){l2(pooKo{p) + lo{<l>o2Ko{p) + 2<j)wKi{p) + k2{p)4>m + <l^2oKo{p)) 



Y [p + im + icloKQ(jp)(l)QQY J 

-^j. ( icloKijp) \ 

010 = jTrp^{f,c} 1 fy ( \A ^^^^^ 

0n = '-Tr, .J { '^^hM^ \ + '_Tr, { ^Mp)K.{p){<Pooh ^ <Po.lo)\ ^ ^ ^^^l) 

4 ' ' \p + im + icloKf){p)4iQ j 4 ' \ (p + im + ic ;oA:o(p)0oo)^ / 
since Ai = 0, 0oi = 

A. ( ichKijp) \ 

4 \p + im + iclQKQ{p)(paQ 



c2 IqKi{p){12(I)qoKo{p) + loi(f'02Ko{p) + 20ioA'i(p) + i^2(p)0oo + 02oA'o(p)) 



(p + im + ic/oA'o(p)0oo) 



2 



^ / icloK2{p) \ , . 

4 yp + iTO + ic/oA^o(p)0oo/ 

, 1 /^^c;i^2(p^_\ 1 /^ C^?oj^o(p)i^2(p)(0OO^l+ 001^0) ^ „ 

4 \p+ im + icloKo[p)(po J 4 y (p + i,m + 2c /oA:o(p)0oo)'' / 

since Ai = 0, 0oi = 

4 \ p + zm + ic/oA^o(p)0oo / 



loK2{p){l2(t)ooKo{p) + ?o(0o2A:o(p) + 20ioj^i(p) + J^2(p)0oo + 02oA^o(p))) \ 
Y (p + + zcZoA'o(p)0oo)^ / 

Thus we finish with a definite set of equations sufficient to determine all the constants 0mn, In- 

1. Evaluate the "ordinary" LO vacuum equation (jB14|B17p . Obtain the values Zo, 0oo 

2. Substitute these values into (|B20IB23p and evaluate 0io,02o 

3. From (|B16p evaluate 0o2 

4. Consider (|B19p as a (linear !) equation w.r.t. I2 and get the value I2 

5. Evaluate directly the remaining values 0i2, 022 

Notice that so far we haven't used any particular form of the size distribution p. The only assumption we did is that 
the width is small, 5p^ <^ 

APPENDIX C: STRUCTURE OF THE The easiest way is to start from (gl]), split the field $ 

0M'''i^''>-CORRELATOR 



In this section we are going to prove that due to the 
chiral symmetry correlator has a form 
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as $ = (tU + 4>' and make the chiral rotation of the quark 
fields 

and will be pa- 



where the matrix U was defined in 
rameterized as 

U{x) = exp (m-r) = 1 + mf — + 



(C3) 



The transformation (|C2p has Jacobian equal to one, so we 
have just to evaluate the rotation of the Dirac operator 
in S. This is convenient to introduce a special notation 



Z Xn 



.x^^F{x)^F,, 



I 



(2^ 



-ip- (x — z) 



(C4) 
(C5) 



i.e. for formfactors the lower index corresponds to dif- 
ferentiation in momentum space (not in coordinate !) 
whereas for all the other quantities we will use the lower 
index for differentiation in coordinate space, e.g. 



U^f_c{x) 



_ du(x) du{x) 



dxi^ 



dx,, 



(C6) 



One more agreement, for the sake of brevity sometimes 
we will not write out explicitly the fiavour dependence 
for the quantities u and a^, implying 



UT 

T' 



(C7) 



Now we are going to evaluate the result of rotation. Note 
that actually we are making a double expansion: up to 
the second order in u ~ and up to the second order 
in derivatives of the field u (assuming to be the first 
order). Rotation of the local part is quite trivial, the 
result is 

[/75,l/2t (^p _^ ^ [/75,l/2t = ... ^ p + + M_^7"75 + 

^ ■ (a^ + u ^) X u + im ^1 - iu-fr, - (C8) 

Notice that the term containing cross-product (in 
isospace) doesn't contribute since it is already a highest 
order in a, u which we keep and it should be combined 
with another vector-isovector combined of u, and the 
derivatives. For the nonlocal part rotation is a bit more 
tricky since we have matrices U at three different points: 



iML{x, z)F{x ~ z)U^'^ {z)L{z, y)F(z - y) + v ^ iMW^^^'^^\x)L{x, z)F{x - z)U^'^ {z)L{z, y)F{z - y)U^^^^'^^^ (y) + v' , 

y' = !Mc/75,l/2,t_^j^(p)[;75,l/2r . $"[/75.1/2^^(p)[^75 4/2,t_ 

a 

I 



where meson fluctuations (I>j and are related by the 
chiral rotation, 

^ ^75,l/2,t$'[/75,l/2,t^ (C9) 

for a moment we will consider only the LO result and 
I 



ignore the fluctuations term v. 



First, we are going to make the Taylor expansion at 
the point z. 



[/75,i/2,t(a.) = ( 1 - m(z)75 



u'^{z)\ f . u{z)u^^,{z) 



+ -«w,p(z)75 



(x - z)^ + 



(a; — z)^{x — z)y + 0(0^, u^, au, d^) 



u,f_i{z)uy{z) u{z)u^y{z) 
4 4 



{z - y)^ + 



(CIO) 



(Cll) 



u^f_c{z)ui,{z) u{z)u^_,y{z) 



{z - y)^{z - y)^ + 0{a^,u^, au, d^) 



Using the general expressions for the nonlocal vertices from the previous Section |^ one can show that 
LF{x - z) = F{x - z) - a^{q)Ffj^{x - 2)75 + ^af_,{q)a^{q)Ffj,,y{x ~ z) + 0{a^) 
LF{z -y)^ F{z ~y) + a^{q)F,,{z - y)j5 + l-a,,{q)a,{q)F^,{z ~ y) + 0{a^) 



(C12) 
(C13) 
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^ W^-^'^-\x)LF{x - z)U^^'^'^{z) = F{x -z)- (^a^(q)75 + «,m(^)75 + \u{z) x + u,^,{z))^ F^{x - z) + 

1 i 

-{a^{q) + u,^,{z)){a,y{q) + w^(z))f)^^(a; - z) + -u^^j5F{x - z) + 0{a^,u^,a'^u,au'^,df_,{udf_,u)) (C14) 

z 

Notice that we drop the terms df^iudf^u) which give just a total derivative 

^ C/75,i/2,t(^)ij^(^ _ y)[/'r54/2(y) ^ _ ^ (^a^(g)75 + ".^(2)75 - \u{z) x (a^(q) + w,^(z))) F^z - + 

1 i 

-{a^{q) + u,f^{z)){dy{q) + w^(z))f)_,t.(z - y) + -u^ulbF{z - y) + 0{a^ , , a'^u, au^ , d^{ud^u)) (C15) 
^ C/75,i/2,t(a.)i(a.^ „ z)U^^^{z)L{z,y)F{z - y)U^''^^'^^\y) = ... = (C16) 

F{x - z)i^(z -y)+ ^0^(9)75 + m,m(2:)75 + ^"(2) X (a^(g) + u,^(z))^ {F^{x - z)F(z - ?/) - F(x - z)Ff,{z - y)) + 

\{a^,{q) + M.,,(z))(a4q) + u^(z))(F^^(x - z)F{z - y) + F[x -~ z)Ff,^{z - y) - 2F^,{x ~ z)F^{z ~ y)) + 
o 

^-u^,^{z)-i^{F^,^{x - z)F{z -y) + F{x - z)F^^{z - y)) + C'(a^ , ...) 

I 



Note that 

1. The term x u doesn't contribute due to its 
isovector structure 

2. The term ~ u^^ may only contribute beyond the 
chiral limit: it contains the second derivative, and 
the "free" u (without derivatives) comes only via 
0(mM)-term. 

3. All the remaining structures have a form {dfj,{q) + 

4. Beyond the chiral limit expansion of Tr In yields 



+0(a^u^m2), (C17) 

where the constants Fij differ only beyond chiral 
limit: 

FL - FL = 2 {Fl^ - Fl^) = -«i m (C18) 

For evaluation of the meson loop contribution (term v') 
one should notice that only the second order, w ^, is es- 
sential. Using expansions (jC14IC15p . one can show that 
in the chiral limit the result depends only upon the struc- 
ture flp + tI/j. Beyond the chiral limit 0(m)-correction 
has a form ~ mM.p(ap + m.^) This completes the proof 
that the correlator of two axial currents has a form (jCip . 



[1] 

[2] 

[3] 

[4] 
[5] 

[6] 

[7] 

[8] 

[9] 

[10] 



J. Gasser and H. Leutwyler, Annals Phys. 158 (1984) 
142. 

T. S chafer and E. V. Shurya k, Rev. Mod. Phys. 70 (1998) 
323 [arXiv:hep-ph/9610451 |. 

D. Diakonov, Prog. Part. Nucl. Phys. 51, 173 (2003) 
(arXiv:hep-ph/0212026 . 

E. V. Shuryak, Nucl. Phys. B 203 (1982) 93. 

D. Diakonov and V. Y. Petrov, Nucl. Phys. B 245 (1984) 
259. 

M. C. Chu, J. M. Grandy, S. Huang and J. W. N egele, 
Phys. Rev. D 49 (1994) 6039 |arXiv:hep-lat/9312071 . 
J. W. Negele, Nucl. Phys. Proc.^SuppI. 73^(1999) 92 
[arXiv:hep-lat/9810053 . 

T. DeGrand, Rev. D 64 (2001) 094508 

[arXiv:hep-lat / OIO6OOI" 



A. G. Williams and J. b. Zhang, Nucl. Phys. Proc. Suppl. 



P. Faccioli an d T. A. DeGrand, Phy s. Rev. Lett. 91 
(2003) 182001 iarXiv:hep- ph/0304219' . 
P. O. Bowman, U. M. Heller, D. B. Leinweber, 



[11] 

[12] 
[13] 

[14] 

[15] 

[16] 

[17] 

[18] 
[19] 



128 (2004) 23 [arXiv:hep-lat/0403002 
M. Cristofor etti, P. Faccioli, M . 
J. W. Negele, [arXiv:hep-ph/0605256| 



C. Traini and 



G. 't Hooft, Nucl. Phys. B 72 (1974) 461. 
D. Diakonov and V. Y. Petrov, Nucl. Phys. B 272 (1986) 
457. 

D. Diakonov, V. Polyakov and C. Weiss, Nucl. Phys. B 



461 (1996) 539 [arXiv:hep-ph/9510232| . 

E. N. Nikolov, W. Broniowski, C. V. Christov, 
G. Ripka and K. Goeke , Nucl. Phys. A 608 (1996) 411 
arXiv:hep-ph/9602274'. 

R. S. Plant and M. C. Birse, Nucl. Phys. A 703 (2002) 
717 ar Xiv:hep-ph/0007340l . 

F. Pena, M. C. Nemes, A. H. Blin and B. Hiller, Braz. J. 
Phys. 29 (1999) 469. 

M. Oertel, arXiv:hep-ph/0012224' 

M. M. Musakhanov and F. C. Khanna, Phys. Lett. B 395 



22 



(1997) 298 [arXiv:he p-ph/9610418\ 

E. D. Salvo and M. M. Musakhano v, Eur. Phys. J. C 5, 
501 (1998) 
M. Musak: 



arXiv:hep-ph/9706537l . 



banov, Eur. Phys. J. C 9 (1999) 235 



a rXiv:hep-pli/9810295]. 

H. Leuhvyler, arXiv:hep-ph/06121 12 

H. Leutwyler, arXiv:0706.3138 [hep-ph]. 

P. Langacker, H. Pagels, Phys. Rev. D 8 (1973) 4595. 

V. A. Novikov, M. A. Shifman, A. I. Vainshtein and 

V. I. Zakharov, Nucl. Phys. B 191(1981) 301. 

C. K. Lee and W. A. Bardeen, Nucl. Phys. B 153 (1979) 

210. 

M. Musakhanov, Nucl. Phys. A 699 (2002) 340. 
M. M. Mus akhanov and H. C. Kim , Phys. Lett. B 572 
(2003) 



H. C 



181 
Kim, 



arXiv:hep-ph/0206233l . 



118 



VI. M usakhanov and M. Siddi kov, Phys. Lett. 
B 608 (2005) 95 [ar Xivihep-ph /041 1 181| . 
R. D. Carlitz and D. B. Creamer, Annals Phys. 
(1979) 429. 

G. V. Dunne, J. Hur, C. Lee and H. Min, Phys. Rev. D 
71 (2005) 085019 a rXiv:hep-t h/0502087 . 
G. 't Hooft, Phys. Rev. D 14 (1976) 3432 [Erratum-ibid. 
D 18 (1978) 2199]. 

S. R. Coleman and E. Weinberg, Phys. Rev. D 7(1973) 



[34] 
[35] 



[36] 
[37] 
[38] 
[39] 

[40] 

[41] 
[42] 



R. Jackiw, Phys. Rev. D 9(1974) 1686. 
P. O. Bowman, private communication (M(m))- 
dependence. See also P. O. Bowman, U. M. Heller, 
D. B. Leinweber, M. B. Parappilly, A. G. Williams 
and J. b. Zhang, Phys. Rev. D 71 (2005) 054507 
arXiv:hcp-lat/0501019 . 

H. C. Kim, M. M. Musakhanov and M. Siddikov, Phys. 

Lett. B 633 (2006) 701 arXiv:hep-ph/0508211 . 

V. Dmitrasinovic, H. J. Schulze, R. Tegen and R. H. 

Lemmer, Annals Phys. 238 (1995) 332 

C. Aubin et al. [MILC Collaboration], Phys. Rev. D 70 

(2004) 114501 [aFXiv: hep-lat/0407028 . 

L. Del Debbio, L. Giusti, M. Luscher, R. Petronzio 

and N. Tantalo, JHEP 0702 (2007) 056 

farXiv:hep-lat /0610059] . 

G. Colangelo, J. Gasser and H. Leutwyler, Nucl. Phys. B 
603 (2001) 125 arXiv:hep-ph/0103088'. 
C. Bernard et al, arXiv:hep-lat/0611024] 
Ph. Boucaud "~ei al. [ETM Collaboration], 



|arXiv:hep-lat /0701012 



[43] A. E. Dorokhov, Eur . Phys. J. C 42 (2005) 309 



[arXiv:hep-ph/0505007 



